Induced superconductivity distinguishes chaotic from integrable bilhards 
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Random-matrix theory is used to show that the proximity to a superconductor opens a gap in the 
excitation spectrum of an electron gas confined to a billiard with a chaotic classical dynamics. In 
contrast, a gapless spectrum is obtained for a non-chaotic rectangular billiard, and it is argued that 
this is generi c for integrable sys tems, 
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The quantization of a system with a chaotic classi- 
cal dynamics is the fundamental problem of the field of 
"quantum chaos" It is known that the statistics 

of the energy levels of a two-dimensional confined region 
(a "billiard") is different if the dynamics is chaotic or 
integrable [^|-^: A chaotic billiard has Wigner-Dyson 
statistics, while an integrable billiard has Poisson statis- 
tics. The two types of statistics are entirely different as 
far as the level correlations are concerned |6|. However, 
the mean level spacing is essentially the same: Parti- 
cles of mass m in a billiard of area A have density of 
states mA/2TTh^ , regardless of whether their dynamics is 
chaotic or not. 

In the solid state, chaotic billiards have been real- 
ized in semiconductor microstructures known as "quan- 
tum dots" 0. These are confined regions in a two- 
dimensional electron gas, of sufficiently small size that 
the electron motion remains ballistic and phase-coherent 
on long time scales. (Long compared to the mean dwell 
time idweii of an electron in the confined region, which 
itself is much longer than the ergodic time igrg in which 
an electron explores the available phase space.) A tun- 
neling experiment measures the density of states in the 
quantum dot, if its capacitance is large enough that the 
Coulomb blockade can be ignored. As mentioned above, 
this measurement does not distinguish chaotic from inte- 
grable dynamics. 

In this paper we show that the density of states be- 
comes a probe for quantum chaos if the electron gas is 
brought into contact with a superconductor. We first 
consider a chaotic billiard. Using random-matrix theory, 
we compute the density of states p{E) near the Fermi 
level {E = 0), and find that the coupling to a supercon- 
ductor via a tunnel barrier induces an energy gap -Egap 
of the order of the Thouless energy Et ~ ?i/ tdwoii • More 
precisely. 



Egap - cNT6/27r, 



(1) 



where N is the number of transverse modes in the barrier, 
r is the tunnel probability per mode, 26 is the mean 
level spacing of the isolated billiard, and c is a number 
which is weakly dependent on F (c decreases from 1 to 



0.6 as F increases from to 1). Eq. (P requires 1 ^ 
A^F <C A/(5, where A is the energy gap in the bulk of the 
superconductor. In this limit p{E) vanishes identically 
for E < -Egap- In contrast, for a rectangular billiard we 
do not find an energy gap in which p — 0, but instead find 
that the density of states vanishes linearly with energy 
near the Fermi level. We present a general argument that 
in an integrable billiard p has a power-law dependence on 
E for small E. 

The system considered is shown schematically in the 
inset of Fig. A confined region in a normal metal (N) 
is connected to a superconductor (S) by a narrow lead 
containing a tunnel barrier. The lead supports N propa- 
gating modes at the Fermi energy. Each mode may have 
a different tunnel probability F„, but later on we will take 
all F„'s equal to F for simplicity. The proximity effects 
considered here require time-reversal symmetry, so we 
assume zero magnetic field. (The case of broken time- 
reversal symmetry has been studied previously [§-0.) 
The quasi-particle excitation spectrum of the system is 
discrete for energies below A. We are interested in the 
low-lying part of the spectrum, consisting of (positive) 
excitation energies £"„ ^ A. We assume that the Thou- 
less energy Et = NY5/2tt is also much smaller than A 

There are two methods to compute the spectrum in 
the regime E, Et ^ A. The first method is a scattering 
approach, which leads to the determinant equation p3| 



Det[l + S'o(£;)S'o*(-£;)] = 0. 



(2) 



The Nx N unitary matrix 5*0 {E) is the scattering matrix 
of the quantum dot plus tunnel barrier at an energy E 
above the Fermi level. Eq. (^) is a convenient starting 
point for the case that the quantum dot is an integrable 
billiard. For the chaotic case, we will use an alternative 
— but equivalent — determinant equation involving an 
effective Hamiltonian [nOl, 



Ho 



(3) 



The M X M Hermitian matrix Hq is the Hamiltonian 
of the isolated quantum dot. (The finite dimension M 



1 



is taken to infinity later on.) Because of time-reversal 
symmetry, Hq = Hq. The AI x N coupling matrix W 
has elements 



(^)'''(2r--i + 2r-v/r3?; 



m = l,2,...Af, n=l,2, 



.N. 



1/2 



(4) 



The energy 5 is one half the mean level spacing of Hq, 
which equals the mean level spacing oiHiiW — 0. 

We now proceed to compute the density of states. We 
first consider the case of a chaotic billiard. The Hamil- 
tonian Hq then has the distribution of the Gaussian or- 
thogonal ensemble ||^, 

P(iJo) oc exp Tri?^) , A = 2M6/tt. (5) 

We seek the density of states 

p{E) = -'K-HuiTv{{E + - H)-'^), (6) 

where (• • •) denotes an average over H for fixed W and 
Hq distributed according to Eq. (||) . The method we use 
to evaluate this average is a perturbation expansion in 
1/M, adapted from Refs. [ p"3p4| ]. Because of the block 
structure of H [see Eq. (||)], the Green function Q{z) = 
{{z — H)^^) consists of four M x M blocks tJn, t/i2, ^21, 
Q22- By taking the trace of each block separately, one 
arrives at a 2 x 2 matrix Green function 



G = 



M 



Trt/ii TtGu 
TrC/21 TrtJ22 



(7) 



(We have multiplied by A/A/ = 26/ n for later conve- 
nience.) One more trace yields the density of states. 



piE) 



To leading order in 1/M, the matrix G satisfies 



G = -y 



nw. 



Z — AGll TTW^ 
2 ^ AG2I Z - 



+ AG12 
AG22 



(8) 



(9) 



where we have abbreviated wi = {WW^)nn- Eq. (|) 
is a matrix-generalization of Pastur's equation |l5). A 
unique solution is obtained by demanding that G goes to 
\/ z times the unit matrix as \z\ — > cxd. 

We now restrict ourselves to identical tunnel probabil- 
ities, r„ = r. For Af > iV > 1/r Eq. (g) simphfies to 



A^GiiJ = 7rzGi2(-Gi2 + 1 - 2/r), 



G2 



Gil, G21 — G12, G]^2 — 1 + '-'II- 



(10a) 
(10b) 

. J. The 

result is that p{E) = for < -Bgap, where -Egap is 



This set of equations can be solved analytically flf 
result is that / 
determined by 
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FIG. 1. Density of states of a chaotic billiard coupled to 
a superconductor (inset), for various coupling strengths. The 
energy is in units of the Thouless energy E'^ — NTS /2n. The 
solid curves are computed from Eqs. (g|) and (p^, for F = 1, 
0.5, 0.25, 0.1. The dashed curve is the asymptotic result ( p^ ) 
for F ^ 1. The data points are a numerical solution of Eq. 
(^), averaged over 10^ matrices Ho in the Gaussian orthog- 
onal ensemble {M = 400, A'' = 80). The deviation from the 
analytical curves is mainly due to the finite dimensionality M 
of Ho in the numerics. 



3k* - 20k^ + 16 



(l-fc)6 (1-fc)^ 

X = Eg,p/ET, k^l- 2/r. 



3P + 8 
(l-k)^' 



(11) 



The solution of this gap equation is the result (|l|) an- 
nounced in the introduction. The complete analytical 
solution of Eq. (0) is omitted here for lack of space. In 
Fig. |l| we plot the resulting density of states. In the limit 
F = 1 of ideal coupling it is given by 



p{E)^ 
Q±{x) 



E^^/?, 
QE5 



[Q+{E/Et)-Q^{E/Et)], 



8-36a;2±3a;V3a!^ 



132x2 - 48 



1/3 



E>Eg,p = 2ET-f 



5/2 



0.6 £;t. 



(12a) 

(12b) 
(12c) 



where 7 = ^(VS — 1) is the golden number. In the op- 
posite limit F ^ 1 of weak coupling we find 

PiE) = E6-\E^ - E^)-'^^, E > Sgap = Et. (13) 



To check the validity of the perturbation theory, we 
have computed p{E) numerically from Eq. (|^) by gener- 
ating a large number of random matrices Hq in the Gaus- 
sian orthogonal ensemble. The numerical results (data 



2 




FIG. 2. Histogram: density of states for a rectangular bil- 
liard (shown to scale in the upper left inset), calculated nu- 
merically from Eq. (^. Dashed curve: Bohr-Sommerfeld ap- 
proximation (|l^). The lower right inset shows the integrated 
density of states, which is the quantity following directly from 
the numerical computation. The energy Et = N5/2tt, with 
N = 200 modes in the lead to the superconductor. 



points in Fig. |1|) are consistent with Eq. (10), given the 
finite dimensionality of Hq in the numerics. 

We now turn to a non-chaotic, rectangular billiard. A 
lead perpendicular to one of the sides of the rectangle 
connects it to a superconductor. (The billiard is drawn 
to scale in the upper left inset of Fig. ||) . There is no tun- 
nel barrier in the lead. The scattering matrix So{E) is 
computed by matching wave functions in the rectangle to 
transverse modes in the lead. The density of states then 
follows from Eq. (||). To improve the statistics, we av- 
eraged over 16 rectangles with small differences in shape 
but the same area A (and hence the same 6 = nfi^ /mA). 
The number of modes in the lead (width W) was fixed at 
N — mv-pW/Trh = 200 (where vp is the Fermi velocity). 
In the lower right inset of Fig. we show the integrated 
density of states v{E) = dE'p{E'), which is the quan- 
tity following directly from the numerical computation. 
The density of states p{E) itself is shown in the main 
plot. 

We have also computed the Bohr-Sommcrfeld approx- 
imation to the density of states, 

PBs{E) =N ds P{s) ^ (5 (i; - (n + ^)TrhvF/s) . 

(14) 

Here P(s) is the classical probability that an electron en- 
tering the billiard will exit after a path length s. Eq. ( p^ ) 
is the Bohr-Sommerfeld quantization rule for the classical 



periodic motion with path length 2s and phase increment 
per period of 2Es/hvF — '!T. The periodic motion is the re- 
sult of Andreev reflection at the interface with the super- 
conductor, which causes the electron to retrace its path as 
a hole. The phase increment consists of a part 2Es/hvF 
because of the energy difference 2E between electron and 
hole, plus a phase shift of — tt from two Andreev reflec- 
tions. For s ^ oo we find P{s) 8{A/Wys~^ , which 
implies a linear i5-dependence of the density of states 
near the Fermi-level, 



PBsiE) 



AE 



2E 
ttEtS ' 



E^O. 



(15) 



In Fig. ^ we see that the exact quantum mechanical 
density of states also has (approximately) a linear in- 
dependence near E = 0, but with a smaller slope than 
the semi-classical Bohr-Sommerfeld approximation. 

We argue that the absence of an excitation gap found 
in the rectangular billiard is generic for the whole class of 
integrable billiards. Our argument is based on the Bohr- 
Sommerfeld approximation. It is known [T^Jl^ that an 
integrable billiard has a power-law distribution of path 
lengths, P(s) s~p for s — > oo. Eq. (|l^) then implies a 
power-law density of states, p{E) oc EP~'^ for E ^ 0. 

To conclude, we have shown that the presence of an ex- 
citation gap in a billiard connected to a superconductor 
is a signature of quantum chaos, which is special in two 
respects: It appears in the spectral density rather than 
in a spectral correlator, and it manifests itself on the 
macroscopic energy scale of the Thouless energy rather 
than on the microscopic scale of the level spacing. Both 
these characteristics are favorable for experimental obser- 
vation. Our theoretical results are rigorous for a chaotic 
billiard and for an integrable rectangular billiard. We 
have presented an argument that the results for the rect- 
angle are generic for the whole class of integrable bil- 
liards, based on the semi-classical Bohr-Sommerfeld ap- 
proximation. It remains a challenge to develop a rigorous 
general theory for the integrable case. 
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